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Abstract 

Arguments are presented to show that in the case of entangled sys- 
tems there are certain difficulties in implementing the usual Bohmian 
interpretation of the wave function in a straightforward manner. Spe- 
cific examples are given. 

The three basic prescriptions of standard Bohmian quantum theory [1] 

are: 

(i) take the wave function ip tohe a solution of the Schrodinger equation, 

(ii) impose the guidance condition p = mdii./dt = V S where S is the 
phase of the wave function ip = Rexp{iS/h), and 

(iii) choose the particle distribution P[tg] at some arbitrary time to (the 
initial time) such that P[f(,] = = Rfto]- This is known in the literature 
as the 'quantum equilibrium hypothesis' (QEH). Given the prescriptions (i) 
through (iii), one can prove complete equivalence between this theory and 
standard quantum mechanics by using the continuity equation for to show 
that P[t] = -R|] for all subsequent times. 

While these prescriptions are self-consistent and work for single particle 
and factorizable many-particle systems, it turns out that non-factorizable 



multi-particle systems sometimes present certain difficulties. This is prin- 
cipally because the velocity equations for such systems are non-separable 
and imply constraints on particle positions which make implementation of 
prescrition (iii) problematic. Although no general method is known yet to 
calculate the trajectories for such systems exactly, analogous to the case of 
the Hamilton- Jacobi theory of many-particle systems in classical mechanics 
P], we give two examples of systems in which the constraint can be calcu- 
lated. 

Let us consider an A^-particle entangled system described by the wave 
function 

tPixi, X2, XN, t) = R{xi, X2, XN, t)ef^3{x^,x,,...,x^,t) 

The guidance conditions which define the velocities are 

dXi 1 dS{Xi,X2, ■■.,X]y,t) 



Vi{xi, X2, -^XN.t) 



dt m dxi 

= fi{xi,X2, ...,XN,t) i = l,2,...,N (2) 



The coordinates not separable, and hence the velocity equations are 

not separable for entangled states, for had they been so, the wave function 
([1]) would have been factorizable which by assumption it is not. Since the 
equations cannot be separated, no general method is known to solve them 
exactly to find the trajectories. 

Consider now the example of a two-particle non-factorizable wave func- 
tion 

where Xi and X2 are the coordinates of the two distinguishable particles of 
the same mass m and momentum p, ^J~L is a normalization constant, and 
a and h are real parameters with a} = \ and a ^ h. Different choices 
of the parameters a, h and p correspond to different wave functions. We 
assume box normalization with —{2N + 1)^ < p{xi — X2)/h < {2N + 1)^, 
being a finite but sufficiently large integer. Then L = (2N + l)7r(a^ -|- 
b^)h /p = {2N + l)7ih /p which is a multiple of the de Broglie wavelength of 
the particles, the only natural scale in the theory. This is a solution of the 
Schrodinger equation provided E = p'^/m. To keep the notation simple, we 
shall henceforth drop the suffix [a, b, p] from the wave function. The position 
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probability density is given by 

R'^{xi,X2,t) = \ij{xi,X2,t)\^ = j[l + 2ab cos2p(xi - X2)/h] (4) 

The important points about this probabihty density, as far as this paper is 
concerned, are that (a) it is stationary, (b) it depends only on (xi — X2) and 
(c) (xi — X2) can take all values in the support of the wave function at all 
times. 

The phase of the wave function ([3]) is 

S{xi,X2,t) = ^arctan[( -) tanp (xi — X2)/^] — Et + rjh 

= h arc ta.n 9 — Et + rjh (5) 

where ^ 

9 = { — --^) tanp {xi — X2) /h (6) 

and 1] = {n + 1)7t for {2n + l)7r/2 < arc tan < {2n + l)7r/2 + 71,1] = —rm for 
-(2n + l)7r/2 < arctan6' < -{2n + 1)tx/2 + tt, n = 0, 1, 2, , ... . In order to 
have a continuous and single-valued wave function, we choose = 0, i.e., the 
principal branch of the function 5'(xi,X2,t) which lies between ±7r/2. The 
Bohmian guidance conditions are therefore 



dxi 1 dS{xi, X2, t) 
dt m dxi 



sec —{xi — X2) 



m[l + eYa + b' h 
p {a-b)/{a + b) 



V2 



m cos^p (xi — X2) /h+ [(a — 6)/(a + b)Y?Av?p (xi — X2) / h 
dx2 1 ^^(xi, X2, t) 
dt m dx2 
p {a — b)/{a + b) 

m cos^p (xi — X2) /h+ [{a — b)/{a + b)]'^sm^p (xi — X2)/h 



(7) 



(8) 



The coordinates Xj are non-separable and so are these equations, i.e., the 
velocity of each particle depends not only on its own position but also on the 
position of the other. Note, however, that ([7]) and ([8]) imply 

Vl + ^2 = (9) 
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which shows that the centre-of-mass of the particles is stationary. 
It also follows from ([7]) and ([HD that 

d{xi — X2) 2p {a — b)/{a + b) 

dt m cos^p{xi — X2)/h+ [{a — b)/{a + b)]'^sm'^p {xi — X2)/h 

(10) 

This equation for {xi — X2) is integrable. Using + 6^ = 1, the solution is 

T , h ab . „ / \ 2j9 ^ . 

[xi - X2) + -— ^ ziTT smzp [xi - X2)/n = — t + fj (11) 



2(a^ — 6^) p (a^ — 6^) m 

where f3 is an arbitrary constant of integration. This is a constraint on the 
particle positions, as we will now show. Since the phase S{xi,X2,t) (Eqn. 
(5)) is well defined for xi — X2 = and the wave function ([3]) does not vanish 
at this point, this equality must hold at some time. Let this time be to- Then 
P = —2pto/m and hence 

1 , h ab . ^ / N 2p 

[xi - X2) + --r-2 sm2p (a;i - X2)lri = — [t - to) (12) 



2(a^ — 6^) p (a^ — 6^) m 

It is straightforward to see from this that Xi — X2 = is the only solution 
at t = to provided Aab < 1. Hence, given Aab < 1, every pair of particles 
in the ensemble at t = to must satisfy this constraint. If there are other 
times tQ,to, etc. at which also the constraint (ITT!) holds, every pair in the 
ensemble must meet at such times too provided 4a6 < 1. Therefore, given the 
range of parameters 4a6 < 1, the particle distribution in the Bohmian theory 
cannot be chosen to match the quantum mechanical distribution (jl]) at these 
times - one must avoid these times to invoke QEH. This makes a Bohmian 
interpretation of these wave functions problematic. No such problem arises 
with single-particle and factorizable many-particle wave functions. 

Another example of a constrained system is the following. Consider the 
two-particle wave function 

^(n,r2,t) = -[- + ] e-i^* (13) 

N riAr2B riBr2A 

where is a normalization factor and 



riA = ^Jxl + {yi-ay + zf r2B = \j xl + {y2 + af + zl (14) 
Tie = + {yi + a^ + z\ r2A = ^/a;^ + (z/2 - a)^ + 4 (15) 
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where the first index i (1,2) in Vij denotes the particle and the second index 
j denotes a point-like slit A of co-ordinates (0,a,0) or a point-like slit B of 
co-ordinates (0,-a,0) which are sources of the two spherical waves in the x > 
space. This wave function is normalizable in a finite volume, analogous to the 
plane wave case. This wave function is separately symmetric under reflection 
in the x axis {yi — > —yi) and the interchange of the two particles 1^2. 
The phase 5* of the wave function f[T^ is (in an obvious notation) 

c, fc ^ riBr2A smk{riA + ^25) + ri^ras smk{riB + rgA) 
b = n arctan 



riBr2A cosk{riA + r2B) + riAr2B cosk{riB + r2A) 

N , , 

= n arctan— (16) 

with 

^ = riBr2A smk{riA + ^2b) + '"ia'"2B smk{riB + r2A) (17) 

D = riBr2A cosk{riA + r2B) + riAr2B cosk{riB + r2A) (18) 

The Cartesian components of the Bohmian velocities of the two particles can 
be computed from S using 

1 ^ ^5* driA ^ dS driB ^ 
m driA dxi driB dxi 
1 ^ dS driA ^ dS driB ^ .^g) 
m driA dyi driB dyi 
1 dS driA ^ dS driB . 
m driA dzi driB dzi 
I ^ dS dr2A ^ dS dr2B 
m dr2A 9x2 dr2B 9x2 
1 dS dr2A ^ dS dr2B . 
m dr2A dy2 dr2B dy2 
1 ^ dS dr2A ^ dS dr2B ^ 
m droA dzo droR dzo 





dxi 


1 dS 


~ "dt 


m dxi 


^yi 


_ dyi 


1 dS 


dt 


m dyi 




dzi 


1 dS 




Ht 


m dzi 




dX2 


1 dS 




~ "dt 


m 8x2 




_ dy2 


1 dS 


dt 


m dy2 




dZ2 


1 dS 


VZ2 


~dt 


m dz2 



(25) 
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where 



dS _ _^ ^2^jj2^-iikriBr2ACOsk{riA + r2B) + r2Bsmk{riB + r2A) 



driA D 
N 

- -^{-kriBr2AS^'ak{riA + r2B)+r2BCOsk{riB + r2A))] (26) 

^'5' _ feri I Ar2/n2i-ir ^^ig^2ACOsfc(riA + r2B)+ri^sinfc(riB + r2A) 

dr2B ~ ^ ^ D 
N 

D 



2\ kriBr2Asmk{riA + r2B) +riACOsk{riB + r2A))] (27) 



The expressions for dS/ driB and dS/ dr2A are easily obtained by the replace- 
ments A ^ B m. the above expressions. These show that the differential 
equations for the velocities of the two particles are non-separable. As we 
have seen, this is a general feature of many-particle entangled systems in 
Bohmian theory. 

It is clear from the velocity equations f|T9|) through (|27|) that the equation 
for each particle can be written solely in terms of its own coordinates provided 

riA = r2B and Tib = r2A (28) 

These are therefore 'integrability conditions' for the velocity equations, or 
equivalently, constraints that the trajectories must satisfy at all times. No- 
tice that no assumption has been made about the initial positions of the 
particles in arriving at these conditions -they are independent of initial con- 
ditions. The existence and properties of other trajectories, if they exist, 
remain conjectural. 

These two examples clearly demonstrate that there are multiparticle en- 
tangled wave functions like ([3]) and flTSl) in standard quantum mechanics for 
which there is no straightforward de Broglie-Bohm interpretation. 

Earlier attempts by Ghose [4| and others [5] to show incompatibility 
between standard quantum mechanics and Bohmian theory were criticised 
mainly on the ground that the initial distributions assumed in these papers 
were incompatible with QEH which is an integral part of Bohmian theory. 
For a full account of the controversies this kind of criticism generated, see 
Struyve and De Baere |6] and references therein. 

An experiment was also performed by Brida et al [7J which claimed to 
simulate the kind of constrained system considered in the second example 
given above. Their claim that the observed results were incompatible with 
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Bohmian theory was criticised for the same reason, namely non-comphance 
with QEH, this time by Oriols [S]. 

The analyses of the two constrained systems considered above in this 
paper show that QEH cannot be invoked for wave functions of entangled 
multiparticle systems in general. Wave functions ([3]) and f|T3l) are examples. 
In the first example the initial conditions cannot be chosen to fit the 
quantum mechanical distribution at arbitrary times, and in the second ex- 
ample (fT3!) the particle distribution is incompatible with QEH at all times, 
independent of initial conditions. 
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